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1 Introduction 



The evaluation of multiloop Fcynman diagrams is currently one of the most important problems in 
Quantum Field Theory. Thanks to the development of analytical techniques in high order perturbative 
calculations, it has become possible to compare precision experimental measurements with the theoretical 
models that try to explain them. There are several of these techniques, and some of the best known can 
even be found in textbooks [1, 2]. One that is not particularly used but which nevertheless has advantages 
for the evaluation of some complicated Fcynman diagrams is called Negative Dimension Integration Method 
(NDIM). The basic foundations of NDIM were initially suggested in the work of Halliday and Ricotta [3], 
which using the dimensional regularization prescription (D = 4 — 2e) make an analytical continuation of the 
dimension D into negative values, something which can be done since the Feynman integrals are in fact 
analytic in arbitrary dimension D. In a previous work [4] we proposed that a more appropriate name for this 
technique should be integration by fractional expansion (IBFE), representing better its mathematical and 
physical basis. 

The purpose of the IBFE technique is to transform the Schwinger integral parametric representation of 
a specific Feynman diagram into an equivalent mathematical structure, which contains several summations 
and Kronecker deltas and which we call multiregion expansion (MRE) of the diagram. This name comes from 
the fact that this particular expansion is made around the values zero and infinity simultaneously, and that 
once the expression is summed using the Kronecker deltas, these expansions become explicitly separated. In 
general the number of summations is bigger than the number of deltas that are generated, and in this case 

cr! 

there are several different ways of summing using the Kronecker deltas, in fact C? — — — different 

01 [a — 0)1 

forms, where 5 is the number of Kronecker deltas and a the number of sums. Each one of these forms give 
rises in general to generalized and multivariable hypergeometric functions. 

Several authors have used and formalized this method, with a wide variety of applications to mainly one 
loop [6, 7, 8, 9, 10] and two loop [11, 12, 13, 14, 15] diagrams. In a previous work we presented a optimization 
of this technique [4], which allows to evaluate more complex diagrams and even certain families of L loop 
diagrams. In this work we will show that is possible to improve the efficiency of the IBFE method when 
it is applied in a modular form to a diagram, that is applying it iteratively to subgraphs or modules that 
contain one or more contiguous loops of the same diagram. The procedure consists in replacing the Schwinger 
integral parametric representation of each module by its corresponding MRE, and simultaneously the diagram 
is reduced topologically into simpler graphs. This reduction allows to get the MRE of the complete diagram, 
as a product of functions that we have called n-loop functions (associated to n loop subgraphs), each of 
them being an MRE related to one of the modules that are present in the diagram. Since each module can 
have different masses distributed in different ways in its propagators, there exist several n-loop functions 
depending on the particular configuration of masses in the propagators. 

The optimization that is reached when applying the IBFE method in this way to a diagram consists in 
a reduction of the obtained MRE with respect to the one that is reached when the diagram is considered 
with all its loops simultaneously. In order to explain the technique in detail, we start with the analysis of the 
simplest topology that can constitute a subgraph of a Feynman diagram: the one loop module (subgraph), 
the bubble diagram, which is going to be the basis that will be used to derive 1-loop functions, useful to 
evaluate a variety of Fcynman diagrams. 

The remainder of the paper is organized as follows. In Sec. II, we describe how to easily obtain the 
Schwinger parametric representation associated to a Feynman diagram, whose mathematical structure is 
the starting point for the application of the IBFE technique. Here we present the basic elements of the 
integration method, and a more rigorous extension is given in the Appendix. In section III we briefly review 
the conventional way to deal with this type of diagrams, which is loop by loop. This will provide a direct test 
for the method proposed here. Later, in Sections IV and V we define and deduce the l-loop functions of the 
bubble module, for all the possible mass and momentum configurations, which are a total of eight functions : 
Gk and Gk {k = A, B, C, D). Finally, in Section VI we develop several applications of the proposed modular 
procedure, which will be compared with its equivalent MRE obtained from an application of the IBFE 
technique considering all loops at the same time. 
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2 Mathematical Formalism 



In what follows we will briefly describe the algebraic elements that are needed to understand the tech- 
nique. 

2.1 Schwinger's pEirametric representation 

The IBFE integration technique is applied directly to the Schwinger parametric representation of a 
diagram, and therefore it is important to have a simple algorithm to obtain this representation, which is 
what we will do in this subsection. 

Let us consider a generic topology G, which represents a Feynman diagram in a scalar theory, and 
let us suppose that the graph has: N propagators, each one associated to the masses {mi, toat}; L 
loops, associated to the independent internal momenta {q'l, Ql}); and E independent external momenta 
{Pi,...,Pe}- 

Using the dimensional regularization prescription we can write the momentum space integral expression 
that represents the diagram in £) = 4 — 2e dimensions: 



G 



jD 



jD 



1 



1 



ITT 2 



r§ {Bl 



(1) 



In this expression the Bj symbol represents the momentum of the j propagator, whose dependence in general 
is given as a linear combination of external and internal momenta. We also define {i/i, z/jv} as the set 
of propagator indices or powers, and which arc considered to have arbitrary values. After introducing 
Schwinger's parametrization it is possible to evaluate the momentum integrals as gaussian integrals, and the 
result of this operation is Schwinger's parametric representation of (1), which in the general case is going to 
be given by an expression of the form: 



G 



(-1)- 



n n^o) ■'o 



exp \ Y.^Xjmyj exp ^-^ 

U 2 



(2) 



For simplicity we have introduced the notation (Tx = dx\...dxN IljLi ^^j' ^) where F is defined as: 



ij Pi-Pj- 



(3) 



The fmiction U and the coefficients Cij are L-lineal and (L + l)-lineal homogeneous polynomials (Symanzik 
or KirchhofF polynomials [16] ) respectively , in tlic^ Schwinger parameters. Both U and Cjj correspond to 
determinants related to a matrix that we have called matrix of parameters M [5] . An important characteristic 
of the coefficients Cij is that they are symmetric {Cij = Cji), which is due to the commutativity of the 
internal product between the independent external momenta associated to the graph G. These determinants 
are defined through the following expressions: 



Cij 



U = 



Ll 



M, 



{L+i)l 



Mil • • • Mil 

Mli ■ ■ ■ Mll 
■■■ Mli 



Mll 



MliL+j) 



M. 



Mr 



(L+i)L M(L+i)(L+j) 



L{L+j) 



(4) 



(5) 
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As can be seen, the determinants are associated to submatrices of the matrix of parameters M, whicli is 
symmetric and of dimension {L + E) x [L + E). The actual form of this matrix can be easily obtained when 
(1) is parameterized and the internal products of all the (internal and external) momenta associated to G 
are expanded, arriving at a quadratic form. The coefficients of such an expansion correspond to the matrix 
elements M^j. For a better understanding of this process let us define for convenience the momentum: 

( Qj if i > i > 1 
0, = <^ (6) 
[ pj_L if {L + E)>j> L, 

with which one can build the {L + E)-vector Q = [Qi Q2 ■■■ Q(^l+e)Y- Using this definition the following 
matrix structure is generated in the integral, after the parametrization application and before the loop 
momenta integration: 



/ d-^ exp / n mh-P - S3 E Q^^^.Q. ' (7) 

and from which we can identiiy the symmetric matrix M. 

2.2 Foundations of the integration method IBFE 
2.2.1 Notation and fundamental formulae 

This technique can be introduced directly by considering the integral expression for the Gamma function. 

The idea is to obtain operational rules, which will allow us later on to work with generalized complicated 
structures of this type of integrals, such as Schwinger's parametric representation of a generic Feynman 
diagram (2). Let us analyze the following integral structure: 

00 

1 If 

J0 " T{0) J '^'^ "^^'^ exp(-Aa;), (8) 



where the quantities A and (3 are arbitrary. Expanding the integrand we get: 

00 

^,=^)Y.<^nA- jdxx^^-\ (9) 

n Q 

where we have defined the factor: 

(-1)" 

Now the integral evaluation will not be done in the usual way, but we define the following operational relation: 

CSO 

j dx X^+^-^ = r ^^^^ h+nfi, (11) 


which makes (9) to be an identity: 
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1 1 ^ (-1)" / , 1 V (-1)" .r /'rf.^^ ^^" + ^^ ^. ' 



A/3 r(/3)^r(n + i) 7 r(/j)^r(n + i) V (-i) 

= ^^"(5;3+n,0 ^^^^ 

n 

For convenience we introduce the following notation: 

dx = (a + /3) , (13) 



where the parenthesis (■) has implicit the constraint associated to the Kronecker delta. This identity is the 
first fundamental formula of the integration method IBFE. With the previously defined notation, equation 
(9) can be written in the following way: 

which corresponds to the multiregion expansion (MRE) of the factor . 

On the other hand, starting from the application of equation (8) and later equation (13) to an arbitrary 
multinomial, wc find the second fundamental formula that the method uses. This expresses the fact that 
a multinomial of a terms can be written as a MRE, in such a way that it contains simultaneously all the 
possible expansions with respect to the ratio of the different terms that are present in the multinomial. This 
MRE can be written in the following form: 

[M + ... + = E - E <u..,n. AT (15) 

where the definition of the factor (10) has been generalized to: 

<^„„..,n. = K-K = (-i)"-"""- r(n, + i).'.r(n. + i) - 

2.2.2 General form of a diagram MRE and its solutions 

Once the parametric representation of the diagram (2) has been obtained, the following step is finding its 
MRE, and for this it is necessary to expand the integrand starting with the exponentials if they exist, and then 
expanding all the multinomials that the procedure is generating according to formula (15). This expansion 
process stops when finally there is only one term, which is a product of all the Schwinger parameters. At this 
point all that is left is to replace the integrals according to formula (13) into its equivalent (.). The result 
is the MRE of the Feynman integral G considered in (1). In the case of a general topology characterized 
by M difi^erent mass scales, A'^ propagators, L loops and a minimal quantity P of invariants associated to 
scalar products of the independent external momenta, then the general form of the MRE of G is given by 
the following expression: 



G = (-1)-- E A iQlr T ft h (17) 

ni ,. . ,nCT 



i=i j=p+i ' k=i r(i'fe) 1=1 r(/3j) 



where we can identify the following quantities: 

• a =^ Multiplicity or number of sums that are present in the MRE of diagram G. 

• A kinematical invariant, which is a quadratic form of the independent external momenta. 
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• aj , j3j , 7j =^ Linear combinations of the indexes {ni, .., rio-}, except for jS^, which is the only one that 
has dependence on the dimension D: 

lii = ^+ni + ...+np. (18) 

The coefficients of the indices of the sum {rij} in the Hnear combinations aj and 7^ are (+1), and in 
the case of /3j the indices have coefficients (—1), except for /J^. 

• N Number of propagators or equivalently number of parametric integrations, which the method 
transforms into N Kronecker deltas. 

• K ==> Total number of MREs performed over the integrand of the parametric representation, which 
in turn generates K constraints or equivalently K Kronecker deltas. This only refers to multinomial 
expansions that are present in the parametric integral. 



In order to find the solutions it is necessary to evaluate the sums that are present in (17), using for this 
purpose the existing constraints between the sum indices, represented by the 5 = {N + K) Kronecker deltas. 
There are several ways to do this evaluation, and in fact the number of different ways to evaluate the MRE 
of G using the Kronecker deltas is given by the combinatorial formula: 

^^^ = W^- ^''^ 
Each of these forms of summing will generate as a result a term in the solution of G, which corresponds to 
a generalized hypergeometric function, and whose multiplicity is given by: 

^i={(T-5). (20) 

In general it is not always possible to use the 5 Kronecker deltas to evaluate a similar number of sums, since 
this will depend on the combination of sum indices for which the sum is going to be done. If this happens, 
these cases simply do not generate contribiitions to the solution. 

In simple terms, the idea of the method is to generate finally an expansion that represents the diagram 
G, the MRE of G, characterized by a multiplicity a in combination with 5 Kronecker deltas. From this it is 
possible to get the solution in terms of generalized hypergeometric functions, series of multiplicity ji whose 
arguments usually correspond to ratios between two characteristic energy scales of the diagram, or in the 
more general case include sums of unitary argument. 

3 Conventional reduction of diagrams with massless bubble inser- 
tions 

3.1 Graphical representation of the reduction 

Many diagrams contain bubble type insertions in the propagators, which in the massless case can be directly 
reduced to a propagator and a factor which is a ratio of Gamma functions. This result is straightforward 
and is very useful for an iterative or loop by loop reduction of topologies or subtopologies built in terms of 
bubble type insertions. The fundamental formula for insertions or bubble modules can be deduced directly 
from the solution of the momentum integral for this diagram: 



dP 

where the indices ai and 02 are arbitrary. The explicit solution of this integral is: 
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where the factor g{ai, 02) is given by: 



g{ai,a2) = {-!)' 



r(ai +a2- f )r(f - ai)r(f - 02) 



r(ai)r(a2)r(£)-ai -02) 
Usually equation (22) is represented pictorially in the following way: 



(23) 




ai + 02 - 



D 



-® 



(24) 



a2 



Another graphical formula, which is useful when two or more scalar propagators are in series, is: 



ai 



a2 



(ai + 02) 



= ®- 



(25) 



This last formula is also valid for a massive theory in the case of two equal masses. This two pictorial 
expressions are enough in order to find the solution of this family of diagrams. Finally the operational problem 
consists in literally reducing series and parallel propagators, and the formulae (24) and (25) represent this 
operation. 

3.2 A simple application 

In order to show the usefulness of the previous pictorial formulae (24) and (25), let us consider the following 
radiative correction to a propagator through which it flows a momentum p: 




(26) 



Then, using equation (24) we can reduce the two bubble type subgraphs that appear in the diagram, and we 
get: 



ai+ 82-0/2 



G = 5(01,02) X 5(03,04) X 




(27) 



aq+aA -D/2 



The resulting loop is easily evaluated after using equation (25), and this finally allows us to get the solution 
of this diagram: 



G = 5(01, 02) X 5(03, 04) X 5(01 + 02 + 05 - ^, 03 + 04 + 06 - X ®- 
which is equivalent to the following expression: 



oi + ... + 06 - ^ 



(28) 
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G = g{ai, 02) X g{a3, 04) x 5(0125 - 0346 - 



2 ^ (p2-)ai + ...+a6-^ 



(29) 



To simplify the notation we have defined for the sum of indices the following symbology: 



aijk... = cii + cij + ak + ■■■ (30) 
Using then formula (23) we finally obtain the explicit solution for the diagram: 



G= (-1) 



so r(ai2 - f )r(f - ai)r(f - a^) r(as4 - f )r(f - a3)r(f - 04) 



T{ai)T{a2)T{D - 012) T{a3)T{a4)T{D - 034) 

(31) 

r(ai23456 - ^)T{D - ai25)T{D - 0346) 1 

r(ai25 - f )r(a346 - T)r(2£' - 0123456) {p'^y^+-+'^e-^ ■ 

4 IBFE and modular reduction of diagrams 
4.1 The bubble module 

The application of IBFE to the bubble module will allow us to find the MREs that generate the 1-loop 
functions which identify the different masses and momenta configurations that can be associated to this 
diagram. We start our analysis by considering this diagram to be associated to the kinematical variable p^, 
and then from this wc can define the following associated 1-loop functions: Gk and Gk (k = A, B,G, D), 
according to the value of this variable and the corresponding mass configuration for this module: 

• Case p^ ^ 

Ga =^ Masslcss propagators. 
Gb One massive propagator. 

Gc => Equal mass propagators. 
Gd Different mass propagators. 

• Case = : This could be obtained from the previous functions, just by taking the kinematical 
variable to zero. Nevertheless, wc; will take a different approach, which uses the series multiregion 
representation obtained directly from the parametric representation of a vacuum bubble. The reason 
is simple, in general the mathematical expression for the MRE is simpler than the corresponding 
expression obtained by taking = in the functions Gk- These functions are useful when diagrams 
associated with vacuum fluctuations are evaluated. The 1-loop functions that we will define for this 
case are analogous to the previous ones: 

Ga => Massless propagators. 
Gb One massive propagator. 

Gc ==^ Equal mass propagators. 
Gd Different mass propagators. 



For our analysis let us start by writing the general integral representation of this diagram in momentum 
space: 

G=[^ I (32) 

J i7r^/2(g2_M2)«.((p-g)2-M2)«'=- ^ ' 
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We shall consider general cases, so the powers aj and are arbitrary. We then obtain Schwinger's parametric 
representation: 



G 



(-1)- 



r(a,)r(afc) 



/ 



ex]){xjMj) exp [xkM^) exp 



dx 



p'' 

Xj + Xk 



{Xj + Xk)'- 



(33) 



Depending on the specific values of the masses {Mj, M^} and of {p^}, we can consider eight possible 1-loop 
functions, which will be discussed in what follows. 

4.2 1-loop functions for ^ 

There are four loop functions that can be defined in this case {Ga, Gb, Gc and Gd), of which only the last 
two are independent . The loop function Ga is a particular case of Gc ov Gd, and in the same way Gb turns 
out to be a particular case of Gd- This will be shown now. 



4.2.1 1-loop function Gi 



Mj ^ Mfe ^ 



Let us start with the most general case, with propagators with different masses. Then the corresponding 
parametric integral is given through equation (33), and the first step is to fractionally expand the exponential 
that contains the invariant p'^, which gives: 



or equivalently: 



G = J2 GD{aj,ak; n; M/, M^) (p^)" = ^ Goiaj, a^; n; Mf, Ml) 



(p2 



where the 1-loop function has been defined as: 

(-1)-* 



GD{aj,ak;n;Mf,M^) = 



r(n mn.] / exp(ar,M/) exp (xkMl) - 

\ J / 



(34) 



(35) 



(36) 



Remember that the factor 0„ is given by (p^ = 

f2 



i-iy 



-. The following step is finding the MRE of the loop 



r(n + i)' 

function GoiajjCik', n; Mj, M|), and for this purpose we expand the mass exponentials in (36) and then the 
binomial in the denominator of the same formula, which gives: 



1 



r(f + n) 



(37) 



Performing the necessary algebra, the MRE for GD{aj,ak', n; M?, M|) is finally given by: 



GD{aj,ak;n;M^,M^] 



(-1)' 



r(aj)r(afe) 



E ^n,m,,m.,l,,l. {'M])-^ {-Mir 



(38) 



with the constraints being defined by the identities: 



Ai = {^ + n + lj + lk), 
A2 = {ttj + n + mj + Ij) , 
A3 = (flfc + n + m/s + Ik) . 



(39) 



9 



The process of replacing the graph in terms of this equivalent 1-loop function (see equation (35) ) can be 
written symbolically through the following graphical formulation: 




(40) 



Derivation of Ga from Go (Case Mj = = 0) : The loop function Ga turns out to be a particular 
case of the 1-loop function Gd- The function Ga can be defined as: 



Ga {aj,ak;n) = GD{aj,ak;n;0,0), 
and whose equivalent MRE corresponds to the expression: 



Ga {aj,ak]n) 



(-1)- 



T{aj)r{ak) 



Ijylk 



A1A2A3 

r(f + n) 



(41) 



(42) 



where the masses were put to zero in (38) and the sums associated to the indices {mj,mk} were eliminated, 
taking later on these indices to zero. The constraints are now given by the identities: 



Ai = +n + lj+lk) : 
A2 = {ttj + n + Ij) , 
A3 = (afc + n + Ik) . 



(43) 



As before it is possible to represent the fractional expansion of the graph symbolically through the following 
graphical equation: 




^ Ga {aj,ak;n) 



-® 



(44) 



ak 



Derivation of Gb (Ccise Mj ^ 0,Mk = 0) : The loop function Gb can be written in terms of the loop 
function Gjj as: 



GB{aj,ak;n;Mj) = GD{aj,ak;n;Mj,0), 
the corresponding MRE is obtained just as for Ga- 



GB{aj,ak;n;Mf) 



(-1)- 



T{aj)r{ak) 



V , , i-M^)^^ A1A2A3 

rrij ,lj ,lk V 2 / 



(45) 



(46) 



The constraints now are given by: 
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A2 = {aj + n + mj + Ij) , 
A3 (flfc + n + Ik) ■ 

The graphical equation which represents the loop reduction is now: 



(47) 




(48) 



4.2.2 1-loop function Gc =^ Mj = Mk = M 

The loop function Gc does not correspond to a particularization of Go- To see this we start with the 
parametric representation already factorized: 



G 



(-1)- 



r{aj)T{ak) 



00 



exp (C/M2) exp (-^p^^ 



U 2 



(49) 



with U = {xj +Xk)- The reason that this loop function cannot be derived from Gd is related to the 
dcnoniinator structure, where the factorization of the polynomial U produces an MRE of (49) which is more 
reduced than when we take Mj = Mk in Gd- 

In equation (49) we expand the exponential which contains the momentum, and obtain the following 
series: 



G = Y^Gc {aj,ak;n; M^) 

n 

where we have defined the 1-loop function Gc {aj, ak\ n; M^): 



(50) 



Gc {aj,a.k;ri; M^) 



(-1)- 



T{aj)T{ak) 



00 

J d'x exp (^{xj + Xk) M"^ 



„n n 



{Xj+XkY 

and expanding the exponential which contains the mass term we have the series: 



(51) 



Gc {aj,ak;n; M^) 



(-1)' 



T{aj)r{ak) 



•^j-^k 



{Xj + Xk) 



¥-+n-li 



Now we shall find the MRE for the binomial in the integrand denominator: 



{Xj + Xk) 



S.+n-U ~ 23 ' 



Jjjk \ 2 



{^+n-li+lj + lk) 



n§+n-k) 



which replaced in (52) and after some algebra allows finally to get the required MRE: 



(52) 



(53) 



Gc {aj,ak;n;M'^) = 



(-1)- 



naj)T{ak) 



li,lj,lk 



A1A2A3 



where the constraints {Ak} are given by: 



(54) 
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Ai = + n h + Ij + Ik) , 
A2 = {uj + n + Ij) , 
A3 = (afc + n + Ik) ■ 

Then the loop reduction in this case is represented graphically by the following pictorial equation: 



(55) 




Notice that the 1-loop function Ga is also a particular case of Gc'- 



(56) 



Ga {aj,ak;n) = Gc {aj,ak;n; 0) . 



(57) 



In all the previous representations or reductions of the bubble diagram one can sec the similitude with the 
conventional equation (24). Nevertheless, the differences are important since in the integration by fractional 
expansion technique (IBFE) massive graphs can be reduced with the same degree of difficulty as in the 
massless cases. 



5 MREs for one loop vacuum diagrams 

The following l-loop functions arc useful when we wish to evaluate certain topologies associated with 
vacuum fluctuations. These functions will be identified as Gk {K = A, B,C, D), in analogy with the ones 
obtained in the previous section for ^ 0. The corresponding MREs of these new cases will be obtained 
directly from their momentum integrals, which corresponds to a bubble without external lines. The resulting 
MRE for each 1-loop function is in fact simpler than the one obtained by making = in the respective 
loop function Gk, that is, it has a more reduced composition of sums and Kronecker deltas. 

5.1 MRE of a bubble with one propagator 
5.1.1 Massless propagator 

This case does not correspond to any of the bubble diagram formulations previously obtained because it is a 
bubble with just one propagator. The momentum integral is: 

a 

[ d^^_ (58) 

J iTT* (9T' 

The representation in terms of MRE is for this case easy to find, using: 

/ ^ D / 2^a = -pT-T / dx x"'^ [ ^ exp(-V) 

J iTT^ (r) r(i^) J J iTT^ 

(59) 

^ ' r(a) ' 




12 



where we have used the Minskowski space integral identity: 



J ITT 2 



exp(-/3g2) 



(-1)- 



(60) 



Finally the required MRE can be written in terms of a pictorial equation as: 




(61) 



One of the properties of Feynman integrals is that they are invariant under momenta scaling, which allows to 
show that the integral in equation (58) vanishes, nevertheless, an MRE has been obtained for this diagram. 
The explanation is quite simple. In fact, first we need to know whether the module constitutes by itself a 
diagram, or it is s product of a reduction of a more complex diagram. For the last option the MRE does 
not vanish, but has only validity within the integration method that is used. This allows to formulate a 
generalization associated to what has been said previously in terms of a theorem: 



Any MRE obtained from a Feynman diagram considering all its loops, 
with a number of constraints (Kronecker deltas) 
bigger than the number of summations, vanishes identically. 



Let us consider a simple example, a two loop vacuum fluctuation in a massless theory, whose conventional 
evaluation is rather obvious, but which nevertheless can be an illustrative example of the theorem: 



G 



The parametric representation is given by: 

(-1)-^ 




G 



d X 



r(ai)r(a2)r(a3) J (^^2-2 + xxxz + X2Xj.y 



Expanding the integrand denominator: 
1 



E 



^ni, 712, 713 



ni +n2 „m +"3 ™"2 +"3 



x\ 



(f + ni +n2+n3) 

r(f) 



(62) 



(63) 



(a;iX2 + a;ia:3 + a:2a;3) ' ni,n2,n3 "^2^ 
and replacing the integral signs by their respective constraints we obtain the following associated MRE: 



(64) 



G = 



(-1) 



-D 



r(ai)r(a2)r(a3) E ' 



( Y + ni + n2 + ns) (ai + rii + 712) (02 + ni + ris) {a^ + n2 + ns) 



r(#) 



(65) 
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The number of constraints (Kronecker deltas) is bigger than the number of summations, and therefore the 
previously mentioned theorem says that it should vanish: 



G = 0. 



(66) 



5.1.2 Massive propagator 

The next degree of difficulty comes from the addition of a mass to the propagator, for which we will use 
two different ways of finding the corresponding MRE. In the first one the usual formalism will be employed, 
which consists in parametrizing and then replacing the resulting integral for its corresponding MRE. The 
second relies in expressing the propagator in terms of its own MRE, and then parametrize the loop integral 
which now does not contain any mass term. 

Alternative I : Direct parametrization of the integral We want to find the MRE of the following 
Feynman integral: 



a,m 




(67) 



Performing the corresponding algebra allows us to obtain the required MRE: 

oo 

I —s-T^ 2^ = FTT / 'exp{xm^) / —5- exp(-a;g^) 

J iiT^ (r - m^) r(a) J J 



(68) 



(-1)- 



r(a) 



Pictorially this is summarized in the following way: 




(69) 



Alternative II : Using the MRE of the massive propagator An alternative form for finding the 

MRE of this module can be implemented using the MRE of the denominator in the loop integral, which 
extracts the mass out of the integral and leaves the equivalent of a vacuum massless bubble. Meanwhile, the 
corresponding MRE of a massive propagator can be easily found using the fundamental equation (15), that 
is: 



(g^ — m^) 



;^Z;2^^ = H '^ni,"2 {-my {qy 



(a + ni + n2) 



(70) 



We can now rewrite the following equation for a massive bubble: 
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</'ni,n. (-m'j X 




(71) 



where the massless bubble can be evaluated using equation (61), and then we finally obtain the MRE for this 
case: 




(-1) ° V- , r_^2\"i (a + ?^i + ri2)(-n2--g) 
r(a) 2- "^"^'^^ ^ ^ r(-n2) 



(72) 



Although (69) and (72) are equivalent, it is clear that (69) is a more compact result, which makes actual 
evaluation easier. 



5.2 MRE of a bubble with two propagators 

5.2.1 Trivial cases : Gc 

There are two cases of vacuum bubbles with two propagators which are reducible to one propagator. This 
happens when the two propagators which form the bubble have equal mass parameters, in which case it is 
possible to use the pictorial equation (25) in order to simplify the graph. 

The first case corresponds to a composition of two massless propagators, which results in an expression 
for the 1-loop function Ga and which is reduced as follows: 



(ai + 02) 



GAiai,a2) 



With the help of equation (61) one then obtains the following expression for Ga- 

ai 





Ga (01,02) = 




(-1) 



.D (oi + a2 — y) 
r(ai + 02) 



(73) 



(74) 



Analogously, for the case of two propagators with the same arbitrary mass m, the resulting 1-loop function 

describes Gc- 
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ai,m 



(ai + a2), m 




Gc (ai,a2;m^ 



02, m 

and then using (69) we get the MRE for this module: 

ai, m 




(75) 



Gc (ai,a2;m^) 




(-1)- 



r(ai + 02) 



)„ (-m j ^ ai + 02 - — + n 



(76) 



02, m 



5.2.2 One massive propagator : Gb 



This case defines the loop function Gb and just as was done before we will consider two different possible 
forms to represent the MRE. This module is given by the integral representation: 



GB(ai,02;TO^) = 




-! 



(77) 



02, m 



Alternative I : Direct parametrization of the integral The parametrization and subsequent evalua- 
tion of the loop integral gives us: 



r(oi)r(a2) 



fd'^q r 
/ exp [- 

J ITT ^ 



ix + y)q^ 



(-1)- 



r(ai)r(a2) 



00 

J dxdy x'^^~^i 



_^ exp(y TO^) 
(x + ?/)"? 



(78) 



After some algebra we get the MRE of this case: 





ai 




GB{ai,a2;m'^) = ® 


a2, m 


J r(ai) ^ ?»1>»2.«3 I ^" j 



(79) 
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where: 



Ai = (f + n2+n3), 

A2 = (ai + n2) , (80) 
A3 = (02 + ni + na) . 

Alternative II : Using the MRE of the massive propagator : The second alternative implies to 
extract the mass of the loop integral, that is, to find the MRE of the massive propagator. Then we have: 



GB{ai,a2;m'^) 




(ai - 712) 



^ ^ ^ni,n2 
m ,712 



11 (02 + ni + 712) 




(81) 



02, m 



or equivalently: 



Ol 

(-l)'"? -e^ /_ 2\"i (Q2 + T^l + ^2) (ai -712- y) 

-Y{^ '^"1."^ i ™ ^ r(ai-n2) 
a2,m 

(82) 

Is evident that equation (82) is the most compact for function Gb- 
5.2.3 Propagators with different masses : Gd 

The next function is associated to a bubble with two propagators of different mass, which defines the loop 
function Gd- Again we will deduce this function in two alternative ways, as was done before: 

Alternative I : Expansion of the complete integral Applying Schwinger's parametrization allows to 
find the following structure for this module: 



GB{ai,a2;m'^) = 




ai,M 




a2,m 

(83) 

Performing the evaluations and replacements already described, we get the MRE for this case: 
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ai,M 




(-1)- 



r(ai)r(a2) „f- 



ra3,n4 



2\"i / ;\/r2\"2 A1A2A3 



r(f) 



a2,m 



where the associated constraints are: 

Ai = (f + 713 + 714) , 
A2 = (ai + 7ii + na) , 

A3 = (a2 + 7l2 + 7x4) • 

Alternative II : Expansion of massive propagators The momentum integral is given by: 

ai,M 



GD{ai,a2;M\m^) 



then we expand each massive propagator: 



ai,M 




1 



a2,m 




a2,m 



(-712 - 714) 



2\ni \ai + 7li + 712) 



r(ai) 



E 



(-"^ ) 



(a2 + 7t3 + 714) 

r(i.2) 




andDIDO using formula (61) we obtain the loop function Gd- 
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ai,M 



GD{ai,a2\M^,w?) = 




02, m 



(-1)- 



ri3,rt4 



A1A2A3 

r(-n2 - n4) 



(88) 



with the foUowing constraints: 



Ai = (ai + ni + 712) , 

A2 = (a2 + 713 + n^) , (89) 

A3 = (-n2 -714- §) . 

In all the massive cases we can use two equivalent ways of finding the MRE of the present module, but 
one of them will produce a mathematical expansion in a more reduced or minimal form. More compact 
MREs allow a considerable reduction of the possibilities of evaluating the sums of the expansion with the 
available Kronecker deltas. Remember that in the case of the MRE of a certain Feynman diagram, each 
alternative way of doing these evaluations generates a term (hypergeometric scries) of the final solution, and 
therefore a minimal MRE allows to eliminate irrelevant or non existing terms (cases where it is not possible 
mathematically to perform the sum using the Kronecker deltas). 



6 Applications 

In this section we will use the previously obtained 1-loop functions and present three examples that will 

show and explain the IBFE methodology applied modularly. Together with this wc will be able to visualize 
the advantages of this procedure compared with IBFE applied to the Schwinger's parametric representation 
of the complete diagram. In our previous work [4] , we show as find explicitly solutions to start of the MRE 
equation of a Feynman diagram. In this work, we just compare the ways to obtain the respective MRE 
associated to a topology. 

6.1 Example I : Radiative correction to the three loop massless propagator 

Let us considcir as a first example the three loop correction to the massless propagator, with momentum p, as 
shown in the graphical formula (90), applying the fractional expansion method loop by loop, and using the 
1-loop functions previously obtained. In this case it is clear that we need only 1-loop functions of type Ga- 



G = ^ 




We start by reducing the two bubble insertions directly, using the pictorial formula (44), which allows us to 
get: 
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(as - ni) 



G= ^ Ga (ai,a2;rii) X Ga (a3,a4;n2) x <8> 




(91) 



(ag - n2) 



We can use (25) to sum the propagator indices, and then apply formula (44) for finally getting the diagram 
MRE in terms of a sum of products of 1-loop functions Ga- 



G= ^ Ga (ai,02;ni) X Ga (a3,a4;n2) X Ga (as - ni,a6 - »^2;?^3) X (g)- 



(92) 



This result can be written also as: 

G= ^ Ga (ai,a2;ni) X Ga (03,04; n2) X Ga (as - J^iifle - "^2; ria) (p^)"^ 

ni,--i"3 

where the 1-loop functions Ga for each case are determined, according to equation (42), by the 
expression: 



Ga (ai,a2;ni) 



(-1)- 



r(ai)r(a2) 



^ ^ 4^ni ,714,715 



A1A2A3 

r(# + ni)' 



with constraints that are given by the identities: 



Ai = ( Y + ni + 'T-4 + n^) , 
A2 = (ai + ni + 714) , 
A3 = (a2 + ni + ns) . 



Similarly one has: 



Ga (a3,a4;n2) 



(-1)- 



E 



A4ASA6 



i^Mr(a4) r(f + n2)' 

^4 = + '^2 + ^6 + nr) , 
As = (as + 712 + "-e) , 
Ae = (a4 + 712 + nr) , 



and finally: 



Ga (as - ni, ae - n2] n^) 



(-1)- 



E 



A7A8A9 



r(as - n.nae - 7^2) r(f + 7^3) ' 



As = (as - 7ii + 7i3 + ns) , 

Ag = (ae - 712 + 713 + »^9) • 

Replacing these series in equation (93) allows to get the diagram MRE: 



G = 



(-1)- 



Ui=i r(ai) „; 



E 



2\n3 



n =1 A, 



T{% + 7ii)r(f + 7i2)r(§ + 713) r(as - 7ii)r(a6 - 712) 



(93) 
following 

(94) 
(95) 

(96) 
(97) 

(98) 
(99) 



(100) 
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In order to show the advantages of modular IBFE, we show in Table I a comparison of the MRE obtained 
in this form with respect to that one that evaluates the complete diagram, that is which uses the parametric 
representation that incudes all loops simultaneously. 



Complete (better factorization) Modular 



Multiplicity multiregion series (cr) 
Kronecker deltas of the expansion (S) 
Multiplicity resulting scries (cr — S) 
Posible contributions to the solution (C^) 



11 9 

11 9 



1 1 



{Table I) 



The expression (100) is equivalent to equation (29), which was deduced for the same diagram but using the 
usual 1-loop function (23). The first conclusion is that the fractional expansion technique is, for cases of loop 
by loop reducible massless diagrams, more tedious than the conventional. In spite of this, we will see that for 
cases in which the diagram contains different mass scales, the present technique es very powerful, and this 
will become clear in the following example. 

6.2 Example II : Propagator with two mass scales and two loops 

Let us consider a graph which consists of a two loop propagator with two different mass scales, as shown in 
the figure: 



^5 

We first reduce the loop associated to the mass m, which uses the 1-loop function Gc (equation (56)), and 
then we get: 



m 



G= ® 




(101) 



(as + 04 - Til) 




(102) 



a5,M 

For the resulting loop reduction we use formula (48), obtaining the massive diagram MRE: 



-n2 




(103) 



ni,n2 



which is equivalent to the algebraic expression: 
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G= ^ Gc (ai,a2;ni;m^) X Gb (as, 03 + 04 -ni; 77-2; M^) 



(104) 



The 1-loop functions Gc and Gb in (104) are given by: 

(-1)-- 



Gc (ai,a2;ni;m^) 
where the constraints are: 



and: 



r(ai)r(a2) 



A1A2A3 



r(f +ni-n3)' 



Ai = + ni - 713 + 714 + "5) , 

A2 = ((Xi + ni + 714) , 
A3 = (a2 + ni + ns) , 



(105) 



(106) 



GB{a5,a3 + 04 - ni;n2;M^) 
with the constraints: 



(-1)- 



r(a5)r(a3 + 04 - 7ii) 



E 



^ ^ r(# + 7i2)' 



"^4 = (-J + "2 + "7 + 7l8) , 

A5 = (as + "2 + ne + 717) , 
Ae = (03 + 04 - ni + n2 + ng) . 

The replacement in (104) allows to finally get the MRE of diagram G: 



G 



(-1)- 



r(ai)r(a2)r(a5) 



E 



(p2)"^ (-m2)"^^(-Af2)"'' A1A2A3A4A5A6 
r(f +ni-n3)r(f + 712) r(a3+a4-7ii) 



(107) 



(108) 



(109) 



Operationally it is not possible to apply directly the conventional l-loop reduction (23) to this massive 
example, as has been done here with the fractionally expanded 1-loop functions. Moreover, it is also possible 
to verify that the diagram MRE obtained loop by loop is far more direct and simple than to get this expansion 
considering all loops simultaneously. This can be seen from Tabic II, where wc observe that the modiilar 
application of the IBFE technique not only reduces the MRE of the diagram in terms of sums and deltas, 
but also brings in a significant reduction of the number of irrelevant terms and therefore the calculation time 
of the solution is optimized. 





Complete (better factorization) 


Modular 


Multiplicity multiregion series (a) 


11 


8 


Kronecker deltas of the expansion (5) 


9 


6 


Multiplicity resulting series {a — 5) 


2 


2 


Possible contributions to the solution (Cg) 


55 


28 



{Tabla II) 



6.3 Example III : Vacuum fluctuation with three mass scales and flve loops 

The next diagram that we will analyze is one without external lines, composed of five loops and three different 
mass scales {mi, 77x2, Af}, arranged as shown here: 



22 



G 




'4 "'2 If "3 



In the first step we reduce the left and right massive loops in terms of the 1-loop functions Gb 

^10 ^11 



G = ^ GB{ai,a2;ni;ml)xGB{a3,a4;n2;ml) 




m m) 



ag 



Applying formula (25) in order to sum the propagator indices, we get the equivalent equation: 
G= ^ GB{ai,a2;ni;rnl) X GB{a3,a4;n2;ml) 





^11+ ^12 -"2 



Then we eliminate the internal bubble: 

G= ^ GB(ai,a2;ni;mi) X 65(03,04; 77-2; mi) X Gc(a7, as; «3;-W^) 





continuing with a reduction of the left loop: 



23 



G= ^ GB(ai,a2;ni;TOi) X G'B(a3,a4;n2;m|) X Gc(a7,a8;'T-3;-M^) 

ni,..,n4 



GA(a9,io - 05,6 - '^3; ^4) x 



-71,4 




ciii + ai2 - n'2 



(112) 



The last loop reduction is represented by the function Ga, which determines finally the MRE of the diagram: 



G = ^ GB{ai,a2;ni\m\) x 05(03, 04; 712; m|) x Gc(a7, og; ns; M^) x 6^(09,10 - «i,«5,6 - «3;«4) 



X GA{aiis2 - n2, -^4). 
The corresponding 1-loop functions in this example are given by: 



(113) 



Gsiaua2;n,,ml) = J.^^^ ^^^^^^^^^^^ (114) 

Al = (Ol +711+715+ 717) , 

A2 = (a2+77l+776), (115) 

A3 = ( Y + J^i + + ne) , 

A4 = (03 + 7l2 + 7l8 + nio) , 

A5 = (a4 + 712 + ng) , (117) 
^6 = {§ +n2+n8 + ng) , 

Gciar, as; n,, M^) = ^A-i^ ^^^^^ ^„3,„^^,..,„^3 j;^^^^^^ A7A3A9, (118) 
A7 = (07 + n3 + nii) , 

As = (08 + 7i3 + 7212) , (119) 
Ag = + "3 - niz + nil + ^12) , 

{ — A10A11A12 

Ga («M0 - ni, as.e - ^3; 774) = r{a, + aio - ni)V{a, + a, - n,) J^^^^ T(fT^' ^''^^ 

Aio = ( Y + + ni4 + Tiis) , 

All = (ag + aw - ni + 7^4 + 7ii4) , (121) 
A12 = (as + 06 - «3 + ^4 + nis) , 

Ga (an + ai2 - na, -714) = (-1)""? ^^7 7, (122) 

r(aii + ai2 - 772 - 7^4) 

{ Ai3 = (an + ai2 - n2 - 714 - f ) . (123) 
Finally we obtain the MRE associated to the diagram: 
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r(a7)r(a8) 11^=1 r(a,) 

1 

r(# + ni)r(§ + n2)r(f + ns - ni3)r(§ + n4)r(f ) 

r(a9 + aio - ni)r(a5 + ae - n3)r(aii + ai2 - n2 - 714) 
Table III shows again a comparison of the IBFE apphed modularly or complete to the diagram. 





Complete (better factorization) 


Modular 


Multiplicity multiregion series (a) 


21 


15 


Kronecker deltas of the expansion (S) 


19 


13 


Multiplicity resulting series (cr — 6) 


2 


2 


Possible contributions to the solution (Cg ) 


210 


105 



{Table III) 



7 Comments 

7.1 Other 1-loop functions associated to the bubble module 

The 1-loop functions for diagrams that contain the bubble type insertions previously deduced do not con- 
template all the possible cases. To see this let us try to solve the following topology: 

a? 
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G = 




(125) 



as 

This diagram contains two massive propagators, each of them characterized by mass M, arranged as shown 
in the figure (continuous line), and the rest of the propagators are massless (dashed lines). It can be easily 
evaluated if Schwinger's parametric representation of the whole diagram is integrated. Let us see: 



, -r -P((-3 + -4) M^) exp ( [^p + (a:.+.2)(x3 + .4)bs A 

_ (-1) f V X1X2 + {Xi + X2) {X3 + X4,) + X5 {X3 + X4) J 

— 5 / S. 
TT T^/_ ^ \xiXo + (x^ + Xo) (x-i + Xa) + (x-i + Xa)] ^ 



G 



n r(a,) •'o 



[xiX2 + (a;i + X2) {x3 + Xi) + x^ {xs + Xi)] = 



(126) 



The integral is writtc'n in a factorizcd form and it is readily evaluated with IBFE. After doing the corre- 
sponding expansions we realize that the obtained MRE has ten summations and nine deltas, so the solutions 
will correspond to one variable hypergeometric functions ( gi^g_i), which is something expected given the 
type of topology and the number of energy scales present in the problem. Moreover, we will have at most 



ten possible contributions that distribute themselves in the kinematical regions 



< 1 and 



M2 



< 1. 
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Now wc will try to solve this diagram, but this time applying the 1-loop functions deduced in this chapter. 
The upper loop is easily reduced using the 1-loop function Ga, getting then the following reduced diagram: 



G = '^Ga {ai,ar,n) x 



(8)->— 1 L^-® 



(127) 



35 

Nevertheless, we obtain one loop which is not possible to reduce with the 1-loop functions previously defined. 

The reason is the presence of the mixed massive propagator in the upper branch of the diagram. Thus it 
is necessary to consider a new loop function that we will call Ge, and which is associated to the generic 
topology of the form: 



G= <g> 




(128) 



Let us try to find the MRE for this new loop function. For this we start with the corresponding Schwinger 
parametric representation. 



G 



(-1)- 



Tiai)Tiaj)T{ak) 



/.I 



exp [xiM-^) exp ( --^^ — p^ 

\ ■^i ~\~ Xj ~\~ 
(.X',; + Xj + Xk) ^ 



When expanding the exponential that contains p^ the following series is obtained: 



G = 'S2Ge {ai,aj,ak;n;M'^) — 



and just as before the 1-loop function is defined: 

(-1)- 



Ge {ai,aj,ak;n;M'^) = 



T{ai)r{aj)r{ak) 



exp (xjM^) x'^ {xi + XjY 



(129) 



(130) 



(131) 



Taking into account the complete expansion procedure finally one finds that the multiregion series for this 
new 1-loop function is: 



Ge {ai,aj,ak;n;M'^) = 



(-1)- 



T{ai)T{aj)T{ak) 



^ y '^n,si,..,S5 

Sl,..,S5 



A1...A5 



r(f +n)r(-n-S2) 



(132) 



with the constraints: 



' Ai = (§+n + S2 + S3), 
A2 = {-n - S2 + S4 + S5) , 

A3 = (fti + Sl + Si) , 

A4 = (flj + S5) , 
A5 = (ttfe + n + S3) . 



(133) 
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Graphically this result can be represented as follows: 



G= (& 




—n 



<^ — >^ — ® ='^Ge {ai,aj,ak;n;M'^) x( 



(134) 



This loop contributes with 6E and 56. Coming back now to our initial problem presented in equation (125), 
and using formula (134) we find the MRE for this problem: 



G = (ai,a2; n) x Ge (as + 0,4, -n, 05; I; M^) (p^)' 



(135) 



which can be compared with the equivalent MRE coming from the parametric representation of the complete 
diagram (Table IV): 



Complete (better factorization) Modular 
Multiplicity multiregion series (cr) 10 9 
Kronecker deltas of the expansion (5) 9 8 {Table IV) 

Multiplicity resultanting series {a — 6) 1 1 

Possible contributions to the solution {Cg) 10 9 

The first we should notice is that in this case once again the expression for the MRE obtained in terms of 
modules is more reduced than the similar expression which evaluates the complete diagram. A second aspect 
is that just as the loop function Ge, there are many others which are variations of the bubble module and 
that have not been included here, although as has been done above any other configuration has an MRE 
which can be trivially found. 



7.2 Generic diagrams with bubble type insertions 

We have shown that the IBFE technique is very useful and simple to apply in diagrams that are built by 
successive one loop insertions, which can be either massless or massive. Nevertheless, the modular reduction 
of a diagram using the IBFE technique is not only readily applicable to this type of topologies, but as will 
be discussed later, the previously deduced formulae are also useful for reducing any bubble type subgraph 
that is contained within a generic graph. As an example let us analyze the following diagram: 




(136) 



Independently of the kinematic characteristics of this diagram, our interest is in the propagator correction 
with a bubble type insertion and with one of the propagators with mass m. Using the respective loop function 
in order to reduce the loop we obtain the following graphical equation: 
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\ 





X 



(137) 



n 



/ 



From a topological point of view one can see that there is a reduction, so to find the MRE of this resulting 
one loop diagram is simpler than the original two loop one. Just as in this example, it is possible to use the 

1-loop functions in generic diagrams whose propagators contain bubble type corrections and its variations, 
which simplifies considerably the search of an MRE of certain diagrams and its solutions. 

7.3 Vr-loop modules and Vr-loop functions 

We have shown that it is possible to apply IBFE to loop by loop recursively built subgraphs and find the 

respective MRE in terms of the 1-loop functions defined above. Nevertheless, the flexibility of IBFE goes 
beyond that, and it is possible to generate also n-loop functions or equivalently module functions that consider 
two or more loops simultaneously. Such a module, with n loops, will lead to a growing number of n-loop 
functions when one takes into account the different possibilities of distributing the propagator masses of the 
n-loop module. In general to work modularly is advantageous for reducing diagrams, since it simplifies the 
search of the MRE. Nevertheless, if the idea is to have less algebraic manipulations and a minimal expression 
for the MRE of the topology, then the previous examples indicate that the MRE is going to be minimal when 
the modules under consideration are compositions of only one loop. 

8 Conclusions 

In this work we have described an efficient and easy way to to implement the integration technique 
IBFE, applied to a specific Feynman diagram. The result that has been reached using the concept of modular 
reduction has two important characteristics: first, in general the diagram multircgion expansion of a given 
diagram is shorter than if the MRE is obtained from the parametric representation of the complete diagram; 
second, the modular treatment operationally facilitates obtaining the multiregion expansion. Let us discuss 
this last point, comparing the modular reduction in cases of loop by loop reducible diagrams with l-loop 
insertions in a massless theory which are evaluated using conventional calculation. The difference is that 
with IBFE it is possible to do a calculation with the same degree of simplicity, although for more complex 
bubble type diagrams it is now possible to include masses. 

The modular reduction idea leads in general to a minimization of the MRE of a diagram, which is more 
evident when the modules to reduce are composed of one loop. In fact, the loop by loop reduction optimizes 
the IBFE application. In more general terms and beyond this work, depending on the one loop topology 
the minimal MRE is reached by a systematic reduction going from the topologically simplest to the more 
complex: bubbles. Triangles, Boxes, etc. 

The importance of a minimal MRE is of course tied to a minimization of the time of finding the solutions, 
since the MRE gets the solutions to the parametric integrals evaluating summations with Kronecker deltas 
in all possible combinations, and now the MRE is shorter (less summations and Kronecker deltas, which 
diminish equally), then there are less combinations that need to be taken into account in order to find the 
solution. Basically what happens is that several combinations which do not lead to relevant terms in the 
solution are eliminated, and given that with IBFE the integration is replaced by a linear system evaluation, 
the minimized MRE implies smaller linear systems, and therefore simpler to analyze. 

The integration technique IBFE is a simple method and at the same time powerful, which can be advan- 
tageously compared with other Feynman diagram evaluation, and which does not require great mathematical 
knowledge. Although here we have concentrated in a scalar theory, this method is directly applicable to other 
theories by previous scalarization of the tensorial loop integrals, which produces a sum of scalar integrals. 
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In general the most complicated case of an A'^ propagator tcnsorial integral is precisely the corresponding N 
propagator scalar integral. This is the reason for improving or creating numerical or analytical techniques 
for the evaluation of scalar integrals in perturbative studies in field theory. 
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Summary of formulae for the bubble module 
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= ^GB{aj,ak;n;Mf) x( 
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Y = {o-j, a-k) 


Composition = 


=^ 0S/1(5. 
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=^ 22/2(5. 
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j =GD{aj,ak;M],Ml) 


ak,Mk 
Composition 


=^ 4S/35. 



B Hypergeometric functions 

In the evaluation of L loop Feynman diagrams, all the possible solutions that can be found are presented 
in terms of multivariable generalized hypergeometric series. The actual value of the variables is given in 
general by ratios of two energy scales associated to the graph [4], including also arguments of value '1'. The 
purpose of this appendix is to provide the necessary information about simple and double hypergeometric 
functions [17, 18, 19, 20, 21], which appear in the solution of the loop integrals, and in particular present 
their convergence conditions, which are what finally define the different kinematical regions associated to the 
solution of a specific Feynman diagram. 
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B.l Definition of the generalized hypergeometric function 

The functions associated to one variable series solutions of a Feynman diagram are always expressible in 
terms of generalized hypergeometric functions, in particular the hypergeometric function of order (g, q — 1), 
which is the only one variable series that is solution to the loop integrals. Thsi function is conventionally 
denoted as: 



,F,_,{{a};{b};z)^ 



q^q-l 



{a} 
{b} 



where the factors {a)/, are called Pochhammer symbols and which are defined by: 

For these functions the convergence conditions are: 

• The generalized hypergeometric functions are given by series defined in the convergence region \z\ < 1, 
while for \z\ > 1 they are defined by analytical continuation. 

• If 1^1 < 1 the scries converges absolutely. Since the variable z represents the ratio between two energy 
scales of the topology, what is obtained is one of the limit representations or solutions, an expansion 
around z ^ 0. 

• If ^; = 1, the necessary requirement for the series convergence is that 3?e (w) > 0, where u) is called 
parametric excess and it is given by the equation: 

• For the convergence in ^ = —1 it is sufficient that (w) > — 1. 
B.2 Some identities of the Pochhammer symbols 

The following identities are very useful for building the hypergeometric function starting from the contribu- 
tions that are obtained of the MRE of an arbitrary diagram G. Aside from formula (139), sometimes it is 
useful to use other identities that are needed when there are factors of the type r(a ± n) and r(a ± 2n) in 
the solutions, such as: 

(a) - r(a-n) _ (-1)" 

^'-^ r(a) (l-a)„' ^^^^> 

^ -"(!)„ (I (-) 

B.3 Two VEiriable hypergeometric functions 

Here we describe two variable hypergeometric functions, which correspond to the Kampe de Feriet generalized 

p:r:u p:r:u 

double hypergeometric function F^-*-" and the generalized hypergeometric F'^'^''" . 
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p:r:u 



B.3.1 Funcion Kampe de Feriet Fi-s-v 
This function is defined as: 



pq^s\v ( ai>-->"p ai,...,ar Ci,...,C„ 
PlT--:Pq bi,...,bs di,...,dy 





p:r:u / 




= Jpq:s:v j 



m {b} {d} 



x,y 



(143) 



p r u 

oo n n n (c^Om „ ^ 

where the convergence conditions of the double series exist if the following relation between the indices are 
satisfied: 



p + r ^q + s + 1, 

p + u^q + v + 1, 

and if furthermore, the arguments fulfil the condition: 

1 1 
\x\(P-<!) + < 1 ,if{p>q), 



(144) 
(145) 



max{|a;| , |y|} < 1 



if (p < q). 



p:r:u 

B.3.2 Function " 



This series, which appears frequently in the solutions to Feynman diagrams, has the following definition: 



^^:s:v I 0;i,...,Q;p fti,...,^^. Ci,...,C^ 

/?!,.. .,/3 bi,...,bs di,...,dy 



" V {/?} W W 



a;, 2/ 



E 



p r u 

n (Q;i)n-m n (%)n H (Cj)m „ 

i=i 3=1 j=i a;" 



(146) 



n! m! 



n n ibj)n n 

i=i i=i i=i 

where the convergence conditions of the double series exist if the following relation between the indices are 
satisfied: 



p + r ^q + s + 1, 
q + u^p + v + 1. 



(147) 
(148) 



All the series that have been found with the technique used here fulfil this indexes condition. The determi- 
nation of the convergence region of the variables can be done using the Horns general convergence theory 
[21]. 
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C Mathematical Formalism of the integration by fractional expan- 
sion Method 

C.l Introduction 

In section (2.2) we already introduced some algebraic aspects of the IBFE technique, which are generated 
from the identity associated to the integral parametrization of loops and which is known as Schwinger's 
parametrization: 

oo 

^ = ^1^^^ ^M-Ax). (149) 



It is possible to find an operational equivalence between the integral symbol and a Kronecker delta, given by: 

j dx ^ r (/3) ^pjyr^ h+nfi- (150) 

For simplicity we have eliminated the integral limits, since this identity has only validity in the context of 
the integrand expansion in (149). This expression is crucial for the development of the IBFE method, since 
in the Feynman diagram evaluation the corresponding Schwinger parametric representation is a generalized 
structure of the expression (149). 

C.1.1 Some Properties 

For the study of some properties of (150) it is convenient to use the notation defined in (13) in order to help 
us formalize the mechanism of the IBFE technique. Then let: 

J dx X-''+'''-^ = {Ul + U2) , (151) 

where ui and 1/2 are indices which can take arbitrary values. 



Property 1. Commutativity of indices : We can explicitly write (151) in two possible forms according 
to formula (150): 



{vi + V2) 



r [Vl) — TTITTT— 



1^2 ^l'l+1^2,0j 



(-1) 

p. . r(^i + i) 



(152) 



Starting from (149) one can show the equivalence of both forms in (152), and for this it is enough to expand 
the integrand exponential and replace (•) for each case: 
Let us consider the following integral representation: 



00 

1 1 /" 

-; — = — r / dx x'^^~^ exp (—Ax) . 



which in terms of a series, in the proposed sense (151), turns out to be: 



A''- + V2) 



(153) 



(154) 
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Selecting now {ui + 1/2) = T {ui) ^ ^^-^^'2'^ ^^1+^2,0^ we directly obtain the equality in (153). Analogously 
selecting (i/i + 1/2) = F (1/2) | ^i^^ ^ ^ui+v2,o ^tnd using then the identity: 



V{y-z) ' ' T{l-y) 

with y = ui and z = 2vi, the equality (153) is finally obtained, which shows the equivalence between the two 
ways of writing {vi + 1/2) ■ 

Another way of writing (151), which is useful for the simplification of terms that contain the factor: 

^ ' (156) 



r(TO + 1)' 

where m is an arbitrary index, is the following: 

Tim + 1) 

{vi + vi) = {vi + U2-m + m) = {-m + m) = T{-m) S^^+^^fi. (157) 

(_ljm 

Notice that explicit use of the Kronecker deltas has been made, in order to simplify the parenthesis (•), 
and since this happens in the context of expansions, the Kronecker deltas remain to indicate the constraints 
between the indices i^i and 1^2- 



Property 11. Significance of the Multiregion Expansion MRE : Let us consider the following 
binomial expansion: 

(Al+A2)^^ (158) 

whore the quantities Ai,A2 and i/ can take arbitrary values. In this case there arc two possible regions or 
limits for the expansion: the region where {Ai > A2) and the region where {Ai < A2). These expansions are 
respectively: 

1. Region {Ai > A2) 



(^.+A.)-.A-gj(ek.(-ii) ^ ,159) 

2. Region {Ai < A2) 



(A + A.)-..-|;j|^(-4ij ^ ,100) 

The factor (z/)„ is the Pochhammer symbol and it is given by: 

= (161) 

We have obtained in this manner expansions in the two possible limits separately. It is possible, however, to 
express both results employing a single series which contains simultaneously both regions. In this sense we 

can say that this typo of expansion corresponds to a multiregion scries representation of the binomial. To 
show this, let us express the binomial (158) using the integral representation of the denominator indicated 
in (149). Then we get: 

00 

(Ai + ^2)^" = J dx x^"-^ exp(-a;Ai) exp(-a;A2), (162) 





35 



and the exponentials are expanded separately, obtaining: 



1 . , ( _'\\ni+n2 r 

t \^ \^ ^' A-n-i An2 / „T''+ni+n2-l 



(163) 



Using the identity (151) we get the multiregion binomial expansion: 

+ = f(?^ ? ? r(n;+i)r(n. + i) ^"^^^^ + + ' ^'''^ 

where according to property (152) we can express the parenthesis (•) in three different ways. In any of them 
we will have the same Kronecker delta which eliminates one of the two sums. Let us see: 

i [Tl^ + ni) C>=Fi/+ni+n2,0 

r(Ti^ + n2) ^-pYjtr^ (5=F^+„i+„2,o (165) 

p. ^ ^ r {Tiy + 1) , 
r (712 + ni) (!»zp^+„i+„,,o- 

On the other hand, the number of possible ways of summing (164) using the Kronecker delta can be found in 
general by evaluating the combinatorial C'^eltos which in this case is = 2. Let us see what happens 
when we sum with respect to one particular index: 

1. Sum respect to n2 

Let us use for this case the following equality: 



+ ni + n2) = < 



and then replacing in (164) we obtain: 

(Ai + A2) - ^^^^^ 2^( 1) ^^^^ ^ A, A, , 

or equivalently: 

which gives the expansion associated to the region {Ai < A2), obtained previously in (160). 
2. Sum respect to rii 

Analogously, we now use the identity: 



(166) 



(167) 



(168) 



V / ^ r(ni + 1) ^ 

[Tl' + Til + 712) = r (=FZ^ + n2) 0^^+ni+n2,0, 



and replacing in (164) gives: 



iA,+A2f'' = At'^Y. 



r(ri2 + 1) V Ai 

712=0 ^ ' ^ 

expression that was found in (159), and valid in the region (^1 > A2). 



(169) 



(170) 
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The fundamental idea that has been exposed in the previous demonstration is that using the definition 
(150) it is possible to make a binomial expansion which differs from the conventional on the sense that the 
expansion is now around zero and infinity simultaneously. This can be generalized for multinomials and 
obtain its MRE: 



The number of different expressions that can be extracted starting from equation (171) is given by all 
the possible forms of evaluating some of the sums, using for this purpose the Kronecker delta generated by 
the same expansion, that is C"' = n; possible forms. Generalizing even more, any function expressed of a 
multiregion series employing a sums and 5 Kronecker deltas, has at most: 



possible ways of being evaluated, and each of these expansions corresponds to a series of multiplicity fj, = 



All the resulting series are series representations with respect to the ratios between the terms of the 
multinomial and all of them correspond to multivariable generalizations of the hypergeometric function. 




(171) 



{a -6). 
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